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Abstract 

Recent results from experiments like PAMELA have pointed to ex- 
cesses of in cosmic rays. If interpreted in terms of Dark Matter 
annihilations, they imply the existence of an abundant population 
of in the galactic halo at large. We consider the high energy 
gamma ray fluxes produced by Inverse Compton scattering of in- 
terstellar photons on such e^, and compare them with the avail- 
able data from EGRET and some preliminary data from FERMI. 
We consider different observation regions of the sky and a range 
of DM masses, annihilation channels and DM profiles. We find 
that large portions of the parameter space are excluded, in partic- 
ular for DM masses larger than 1 TeV, for leptonic annihilation 
channels and for benchmark Einasto or NEW profiles. 

1 Introduction 

While compelling evidence for the existence of Dark Matter (DM) now comes from a num- 
ber of astrophysical and cosmological probes (such as galaxy rotation curves, weak lensing 
measurements and the precise data from the Cosmic Microwave Background observations 
and the Large Scale Structure surveys of the Universe [1]), no explicit detection has been 
confirmed yet. The indirect detection strategy relies on the possibility of seeing signals of 
the presence of DM in terms of the final products (e^, p, d, j, u . . .) of DM annihilations 
in the galactic halo, on top of the ordinary cosmic rays. 

The recent positive results from a number of indirect detection experiments have sug- 
gested the possibility that indeed such a signal has been seen. In particular, the signals 
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point to an excess of electrons and positrons. The PAMELA satellite [2] has reported a 
significant excess above the expected smooth astrophysical background and a steep rise 
of the positron fraction e'^ /{e^ + e~) above 10 GeV up to 100 GeV [3], compatibly with 
previous less certain hints from HEAT [4] and AMS-01 [5]. At the same time, no signal 
in the p fiuxes is seen, up to the maximal probed energy of about 100 GeV [6], posing 
stringent constraints on Dark Matter models. Indeed, as discussed in detail in [7], these 
results, if interpreted in terms of DM annihilations, imply either 

(i) a DM particle of any mass (above about 100 GeV) that annihilates only into leptons 
(DM DM e+e-,/i+/i-,r+r-) or 

(ii) a DM particle with a mass around or above a few TeV, that can annihilate into any 
channel (i.e. DM DM W~^W~ , ZZ,bb,tt, light quark pairs and the leptonic chan- 
nels above) possibly producing anti-proton fiuxes at energies above those currently 
probed by PAMELA. 

In any case, a very large annihilation cross section is needed: of the order of lO^^'^cm'^/sec 
up to 10~^° cm^/sec or more, depending on the mass of the candidate and the annihilation 
channel [7]. These numbers are much larger than the typical cross section required by 
DM thermal production in cosmology (~ 3 ■ 10~^^cm^/sec [1]). They can be justified in 
specific models in terms of some enhancement mechanism which is effective today but 
not in the early universe, such as a resonance [7, 8] or Sommerfeld (see [9, 10, 7], and 
then [11, 12]) enhancement, the presence of an astrophysical boost factor due to DM 
substructures (unlikely [13]), or a combination of these (see e.g. [14]). For our purposes, 
they are just an input required by data. 

An excess in the flux of + e~ has also been reported by the ATIC [15] and PPB [16] 
balloon experiments at about 500-800 GeV. The HESS Cerenkov telescope, too, has pub- 
lished data [17] in the range of energy from 600 GeV up to a few TeV, showing a steepening 
of the spectrum which is compatible both with the ATIC peak (which cannot however be 
fully tested) and with a power law with index —3.05 ±0.02 and a cutoff at ~2 TeV. If inter- 
preted in terms of DM annihilations, the ATIC peak would clearly pin down the DM mass 
at about 1 TeV, and as a consequence select the possibility (i) above [7]. These balloon 
and Cerenkov signals will soon be checked by the FERMI satellite that, despite its main 
mission as a gamma ray telescope, will be able to measure e"*" + e~ with unprecedented 
statistics [18]. ^ 

Of course, the origin of these excesses could simply lie in ordinary (albeit possibly 
peculiar) astrophysical sources, such as one or more pulsars [19], sources of CR in galactic 
spiral arms [20], aged SuperNova remnants [21] or exploding stars [22]. In this case, the 
sources would be located in the galactic disk, and moreover not too far from the Earth, 
since quickly loose energy when travelling from more that about 1 kpc away. These 
explanations will be confirmed or ruled out by further, more precise measurements of the 
spectra and possibly improved computations of the expected yields. 

If instead Dark Matter annihilations are at the origin of the excesses, then an abundant 
population of high energy is being created everywhere in the DM galactic halo. In 

^Note added: The FERMI and the HESS collaborations have published new data in [55] and [56] after 
this work was completed. Their inclusion in the DM fits, replacing the ATIC data with which they are 
in tension, does not modifies much, however, the preferred DM green and yellow regions in fig. 3 (see 
e.g. [29]). Indeed: (i) such new data are in agreement with the PAMELA results and (ii) the HESS points 
show indications for a smooth cutoff that still pins down the DM mass at around a few TeV. Our analysis 
here remains therefore essentially unchanged. 
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this case, an inevitable associated signal consists of high energy gamma rays (from tens of 
MeV up to the DM mass, which can reach several TeV) produced by the Inverse Compton 
Scattering (ICS, reviewed below) of the photons of the diffuse radiation fields in the galaxy 
(the CMB, the InfraRed -IR- background and the optical starlight -SL-) on these energetic 

These gamma ray signals have to be compared with the results of the EGRET satellite, 
that published data relative to several windows of observation of the galactic diffuse 7-rays. 
Perhaps more importantly, the FERMI satellite has started publishing new results [23], 
checking in particular the EGRET findings and improving very significantly the statistics. 

The aim of this work is to see whether these existing and upcoming gamma ray ob- 
servations have the power to test or constrain the DM interpretation on the basis of the 
expected associated ICS radiation. 

In terms of their relation to DM annihilations, the ICS signals have been previously 
considered e.g. in [24, 25, 26, 27], in specific cases and specific DM models. Our ap- 
proach differs however on a number of points: we carry out the analysis in a DM-model 
independent way (by considering a set of primary DM annihilation channels, scanning a 
wide range of DM masses and addressing the astrophysical uncertainties, e.g. varying the 
DM galactic distribution profiles), we consider several observational areas/datasets, and 
we perform most of the computations ourselves in a semi-analytical way, as opposed to 
running dedicated codes such as GALPROP [28] . In particular, this last point implies that 

(a) we adopt the simplifying assumption of neglecting the diffusion processes of and 

(b) we adopt simplifying analytical approximations for the density fields of IR and optical 
light, modeled on the detailed results in the literature (see the discussions below). While 
the full numerical approach is of course more precise and probably time effective [29] , our 
custom semi-analytical computation can be quickly adapted to study a number of cases 
and allow us more control on the DM input models. When possible, we have validated our 
results with the results obtained with GALPROP in [24, 26], finding a good agreement. 

The rest of the paper is organized as follows. In Sec. 2 we review the details of the com- 
putation of the ICS signal on the basis of the DM ingredients. In Sec. 3 we briefly describe 
the regions of the sky on which we focus our analysis, the corresponding observational 
datasets and the implementation of the formalism discussed before. In Sec. 4 we present 
and discuss our results and their implications for DM. Fig. 3 shows our main results. In 
Sec. 5 we draw our conclusions. 

2 Basics of Inverse Compton Scattering computations 

The Inverse Compton Scattering process consists of low energy photons that upscatter in 
energy on high energy electrons and positrons (as opposed to the conventional Compton 
Scattering of high energy photons on electrons at rest [30]). A pedagogical review is found 
in ref. [31], on which we will base our discussion and notations. 

In the present case, the high energy are produced by DM annihilations in any given 
point of the galaxy, with a density determined by the DM distribution profile and with a 
spectrum dictated by the primary annihilation products (as discussed below). The bath of 
low energy target photons consists of three main contributions: (i) the starlight originating 
from stars of the galactic disk (at optical wavelengths), (ii) the infrared radiation produced 
by the absorption and subsequent re-emission of such starlight by the galactic dust and 
finally (iii) the ubiquitous microwave photons of the CMB. 
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We will review in the following subsection some basic details of the computation of 
the ICS signal. The uninterested reader finds an executive summary in sec. 2.2 and the 
practical details of our implementation of such formalism in sec. 3. 



2.1 Derivation of the ICS gamma-ray spectrum 

We want to compute the differential flux d^/dei of high energy photons produced by the 
ICS processes, coming from an angular region of the sky denoted by AVt. Here and in the 
following we will denote the energy of the scattered photon with ei and the energy of the 
original photon as e, in the reference system of the photon gas. Such flux is determined as 

f = i/dn/ d.-''^" (1) 

dei ei J Jiinc-of-sight 47r 

in terms of the emissivity j(ei, r) of a cell located at a distance r from the galactic center. 
The coordinate s runs along the line of sight joining the observer at Earth with the point 
r. In general, for any radiative process, the emissivity is obtained by a convolution of the 
spatial density of the emitting medium with the power that it radiates [32]. In our case 
then 

/■a/dm 

j(ei,r) = 2 / dEViei,E,r)n,ir,E), (2) 

where P(ei, E, r) is the differential power emitted into photons of energy ei by an electron 
with energy E and ne(r, E) is the number density at r of electrons with such energy. The 
minimal and maximal energies of the electrons are determined by the electron mass me and 
the mass of the annihilated DM particle Mdm- The overall factor of 2 takes into account 
the fact that, beside the electrons, an equal population of positrons is produced by DM 
annihilations and radiates. 

The quantity ne(r, E), in full generality, has to be determined by solving the diffusion- 
loss equation in each point of the galaxy [33, 34] 



1^ 

j,2 Qj. 



or 



+ v^-l - ^^ir\)/-l + 1^ Wf] = (3) 
dr 3r2 dr dp p"^ dp Anp"^ ' 

V ^ ' '■^^ ^ V ' ^ V ' ^ V ' 

diffusion advoction convection radiative losses source 

written here in terms of / = ne/(47rp^) with p the electron momentum. The advection 
(convection) term describes the flow, of velocity v, towards (away from) the Black Hole 
at the Galactic Center. Since we will be interested in portions of the galactic halo that 
are well outside of the BH accretion region, of typical very small size -Racer — 0.04 pc, 
we assume that these advection and convection terms can be neglected. We neglect the 
diffusion term as well, in order to be able to solve eq. (3) in a semianalytic way (see below). 
This is a more drastic approximation since (for the electron energies that we consider) the 
characteristic time due to radiative losses can be of the order of 10^^ seconds [34], a time 
in which relativistic electrons can diffuse over a distance of about one kpc comparable to 
the size of the smaller regions that we will consider. Diffusion could therefore redistribute 
somewhat the population of electrons and positrons that produce the ICS signal. We will 
verify a posteriori (see Sec. 4) the quantitative impact of this effect and the importance of 
our approximation. Finally, in the radiative loss term we include only the ICS process. 
In principle all processes should be included: synchrotron radiation, bremsstrahlung, ion- 
ization and ICS. However, for relativistic electrons and for the typical magnetic field that 
are present in the galactic halo {B ~ l/xG), ICS turns out to be the dominant one [34]. 
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In short: the emission due to ICS is the dominant process suffered by the electrons and 
positrons produced by DM annihilations, and we assume that they do not diffuse away from 
the production point before having radiated away most of their energy. In this regime, the 
equation is simply solved by 

n,{r,E) = j—- dEQ,{E,r), (4) 

£{E,r) Je 

where S{E, r), related to the p of eq. (3), denotes the total (i.e. into photons of any energy, 
not to be confused with V) rate of electron energy loss due to ICS. 
The source term due to DM annihilations is simply given by 

Q'(E.r) = M^—, (5) 

where {uannv) is the DM annihilation cross section, p(r) is the galactic DM distribution 
profile and dN^/dE is the spectrum of the produced electrons, per annihilation. 

The derivation of S{E,r) and V{ei, E,r) requires some effort, but is straightforward 
in terms of known Compton kinematics [31]. The total energy loss rate for an electron of 
energy E is given by 

in terms of the rate of scatterings dNE,e/ dtdei on photons of energy e into photons of energy 
ei, times the energy lost in a scattering (ei — e) and integrated over all initial and final 
photon energies. The quantity dNE,e/ dtdei is more conveniently computed by moving to 
the frame of the center of mass^ as 

dNEe ff dNEe dt'de\ , , 

' - " ' ^ (11) 



dt dei JJ{n[,e') dt' de[ d^l[ de' dt dei 



^Compton scattering kinematics. In the frame of the center of mass (where the electron is at rest and 
the photons appear to be incident on it) all quantities are denoted with a prime: e' denotes the initial 
energy and e[ the final one. One has 

e' = 7e(l - cos 9), ej^j^ = ^, e„,^^ = 2je, (7) 

where 9 is the approach angle of a photon in the photon gas system and 7 S> 1 is the Lorentz factor. One 
also has the well-known relation 

e' 

= 1 + ^(1 -cos9[y 

where 9[ is the scattering angle in the center of mass. The relation between the final energies in the two 
frames reads 

ei = 7e'i(l-cos^'i). (9) 

The limit of Compton scattering in which the energy of the photon before scattering in the electron rest 
frame is much less than the electron mass defines the Thomson limit. In this case the final energy of the 
photon, in the same frame, is unchanged (the electron looses a negligible fraction of its energy at each 
scattering). Back in the photon gas frame, these conditions translate in the fact that the initial photon 
energy is upscattered to a maximum energy that is larger by a factor 47^ . In formulae 

e'<m, e'i~e', £'"^'^ = 4672 [Thomson limit]. (10) 

It is easy to check that, for typical electrons of DM origin {E < 1 TeV, 7 < 2 10®) the Thomson limit is 
verified for the scattering of CMB photons, since ej^^x = 2 7ecMB ^ 100 eV ^ m, but not for starlight 
photons of typical temperature 0.3 eV, for which e^^.^ — 1 MeV. 
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with dt'/dt = 1/7 and de'i/dei ^ 1/[7(1 — cos^']^)], having denoted the Lorentz factor of the 
electron (always assumed to be highly relativistic) as7 = i^^/?Ti^l. In the new frame, 
the rate of scatterings is given by 

dNE.f dn'(e';e) da ^^^^ 



dt' de[ dn[ de' de' de[ dVt\ 

Here the last term is just the differential Compton cross section, and dn'{e'; e)/de' expresses 
the (differential) density of photons within de' due to photons that are within de in the 
photon gas frame. This quantity, in turn, is to be determined on the basis of the density 
of photons in the original reference frame, which can be done by using proper relativistic 
invariants. One obtains [31] 

1 e' 
dn'{e'; e) de' = -n{e, r) de d{cos 9) — , (13) 

where finally n{e,r) is indeed the photon distribution in the photon gas frame. The last 
formula holds if one neglects possible anisotropics in the spatial distribution of the target 
photons. This is of course precisely true only for CMB photons while it is an approximation 
for starlight and infrared photons. 

The Compton cross section is given, in full generality, by the Klein-Nishina formula 

da 3 /e'^^ 



de'-^dQ'^ 167r 



^t(^1 f^ + ^-sin^^'iVfe'i ^ ), (14) 

' "7 Vi e' 7 l + ^(l-coser' 



where a^ = Snr'^/S = 0.6652 barn is the total Thomson cross section in terms of the 
classical electron radius Vg. In the Thomson limit (of eq.(lO)), it simplifies to 



da 



rfe'i c/fi'i 167r 



(Tt(1 + cos^ ^'1)5 (^e'l - e') [Thomson limit]. (15) 



It is now possible, by solving the integrals in the electron rest frame, to compute the 
rate of scattering of eq.(ll) as [31] 



dNE,e _ o n{e,r)de 
dtdei ~ 472e 

in terms of 



2 



2glng + g + 1 - + \Y^q^'^ ~ 



(16) 



'i=^^J^^V withr, = ^, 61 = ^, in-1^^0<g<l. (17) 
1 e(l — eij m 7m 47"^ 



In the Thomson limit, the same expression holds with ^ (so that the last addend 
vanishes) and, since ei <^ 1, g ^ ?/ = So finally, plugging this result into (6) one 

obtains explicitly'^ 

£{E,r) = 



/•OO /» i 

3(Jt / de e I dq n{e, r 
Jo Jl/472 



:i + r.g)3 



2q\nq + q + l-2q^ + ^^^^(1 - q) 

2 1 + i 



:i8) 



^The integral in q can be done analytically, it just leads to a long expression. 
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In the Thomson hmit it reduces to the particularly compact expression 

£{E,r) = -crT7^ / de e n{e,r) [Thomson limit]. (19) 
3 Jo 

that shows that the total energy lost is proportional to the second power of the electron 
energy E = 7m and to the energy density of the photon bath. 

The computation of V{ei, E,r) proceeds in the same way, apart from the fact that no 
integral over ei is involved (by definition of V). One has in fact 

P(,.E.r) (20) 



from which 

VieuE,r) 



2 1 e? . (21) 
2glng + g + l-2g2 + ^_(l-g) 

2 1 — ti 



Ji/aA Ml^il-ei 
In the Thomson limit 

P{ei,E,r) = —^ei dy ^ ^ [2y\ny + y + 1 - 2y'^] [Thomson limit]. (22) 

47 Jo y 

2.2 Summary 

With the ingredients above, the differential flux of ICS photons (of energy ei) from a region 
AQ of the sky is given by 

^JL ^ L(^r,J^ f ,n [ f^ff^MV r~£^(!^r(£), (23) 

Ai £1 4ir "Ml„]^n Ao... % V Ce / -L S(E, r) 

where Y{E) = J^^^ dE ^ is the number of electrons generated with energy larger than 
E in one annihilation. Here ~ 8.5 kpc is the distance of the Sun from the galactic 
center and ~ 0.3 GeV/cm^ the local normalization of the DM density. This expression 
is subject only to the assumption of highly relativistic electrons and isotropy of the target 
photon gas. 

The total energy loss S{E,r) is given in eq. (18), the power V{ei,E,r) in eq. (21). 
In the Thomson limit (valid e.g. for TeV electrons scattering on the CMB, but not on 
starlight), eq. (19) and (22) can be used. 

The r dependence in V and S comes from that of the photon bath n(e, r). If it is not 
present (such as for the case of the uniformly distributed CMB photons) or it is neglected 
(as we will assume below for starlight and infrared photons), the computation of the flux 
simplifies considerably since the integral over r can be performed separately and the result 
cast in terms of the usual average geometrical factor J [35]. Explicitly 



dei ei 47r M^^ J^ £(E) J An Ji.o.s. "^e \ Pq 



(24) 
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pcri nn 


longitude I 


T 

u 

IsoT NFW Einasto 


A/sL Mr Mcmb 


'5x30' 


0.25° < \h\ < 4.75° 
330.25° < I < 359.75° 
0.25° < I < 29.75° 


10.0 52.1 95.5 


1.7 -10-" 7.0-10-5 1 


1 U A UU 


300.25° < Z < 359.75° 
0.25° < I < 59.75° 




9 7.1 n~i2 7 n . 1 i 


'10-20' 


10° < \b\ < 20° 
0° < / < 360° 


2.3 3.4 4.3 


8.9 -10-" 1.3.10-5 1 




T, = 0.3 eV 3.5 meV 2.725 K 



Table 1: Summary of the observational regions that we consider, together with the corre- 
sponding values of the average J factor for different DM halo profiles and the parameters 
of the modelization of the Inter-Stellar Radiation Field. 

3 Datasets and Application 

A number of observations of the diffuse gamma ray emission are available in the literature. 
In particular, the Energetic Gamma- Ray Experiment Telescope (EGRET) instrument, 
operational in the 90 's aboard the Compton Gamma Ray Observatory, has produced data 
on several windows sitting across the galactic plane [36] or outside of it [37] . The FERMI 
satellite [38] has started exploring the same areas, with much higher sensitivity, and other 
regions on the sky. 

On the basis of the discussion in [39], we consider for definiteness the following three 
regions (see the summary in Table 1): 

• A rectangular region encompassing the inner galaxy area, at galactic latitude 0.25° < 
|6| < 4.75° and galactic longitude 330.25° < / < 359.75° and 0.25° < / < 29.75° 
(denoted in the following as '5x30' region); we take the EGRET data points from 
the reanalysis in [40]. They reach about 80 GeV. 

• A rectangular region larger then the previous one, at galactic latitude 0.25° < |6| < 
9.75° and galactic longitude 300.25° < / < 359.75° and 0.25° < / < 59.75° (denoted 
in the following as '10x60' region). This is the original area considered in [36]. We 
take the EGRET data points from [39] . They also reach about 80 GeV. 

• The intermediate latitude region defined by 10° < |6| < 20° in galactic latitude, at all 
galactic longitudes (denoted in the following as '10—20 strips'); we take the EGRET 
data from [39]"^. The datasets in this region extend up to only ^ 10 GeV in energy. 

For this region we also consider the preliminary data of the FERMI satellite, on the 
basis of [23]. These data do not confirm the presence of the 'GeV excess' reported 
by EGRET, but instead follow the expected purely astrophysical background. Given 
the preliminary nature of the FERMI results and the fact that they are currently 
provided for this region only, we will here show exclusion curves for both datasets. 
But we comment later on the implications of the possibility that FERMI does not 
confirm the 'EGRET excess' anywhere. 



■^The original reference [37] focussed on a larger portion of the sky. 
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Figure 1: Left: The Dark Matter galactic distribution profiles considered in the text. 
Right: Modelization of the total Inter-Stellar Radiation Field (ISRF): the dashed lines 
reproduce the computations of [46, 4V- From top to bottom: Galactic Center (R = 
kpc, z = kpc - black), R = 4 kpc, z = kpc (blue), 5 kpc outside of the galactic 
plane (magenta). The solid lines represent our modelizations in terms of (renormalized) 
blackbody spectra for the StarLight and InfraRed photons. 



In order to compare with these datasets the ICS signals discussed in the previous section 
we need to specify prescriptions for the electron spectrum dNg/dE, the distribution of 
target photons n{e,r) and the distribution of Dark Matter in the galaxy p(r). 

We consider the electron spectra from DM annihilations for the following primary chan- 
nels: 

DMDM ^ e+e-,fi^fi-,T+T-,W^W-,bb and ti, (25) 

computed with the use of the PYTHIA MonteCarlo code [41] as described in detail in [7].^ 
We scan over a large range of DM masses, from 100 GeV up to several tens of TeV. 

As for the photon density distribution n{e), we adopt the following prescriptions, based 
on the precise modelizations of the Inter-Stellar Radiation Field discussed in [46, 47]. First, 
for each region that we study we consider an average density field that does not depend on r 
(but is different for each region, as discussed below). This allows us to use the formalism of 
eq. (24) to compute the ICS signal. Second, we choose to always model the total radiation 
density as a superposition of three blackbody-like spectra (i = 1, 2, 3) 

^-(^)=^^^(^^7^r3T) (26) 

with different temperatures: one for the CMB with Tcmb = 2.753 K (which is of course 
not an approximation), one for the InfraRed light with Tin = 3.5 ■ 10~^ eV and one for 
the StarLight with Tsl = 0.3 eV. Fig. 1 shows that this approximation is actually quite 
good, when compared with the precise computations, in grasping the general shape of the 
spectrum (except perhaps for the details at the junction wavelengths between starlight and 
infrared light, where the density is however suppressed). In order to take into account the 

^We stick to the case of "direct" annihilation of DM particles into a pair of Standard Model particles, 
i.e. we do not study models in which the annihilation proceeds into some new light mediator state (see for 
instance [29, 11, 42, 43, 44, 45]). These "cascade annihilation processes" lead generically to softer electron 
spectra. 
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spatial variation of the SL and IR fields of radiation throughout the galaxy, we normalize 
differently the blackbody spectra in each of the three regions that we study, with the 
overall coefficients Ma {a = 1, 2, 3). More precisely, for the '5x30' region we take the ISRF 
characteristic of the Galactic Center (i? = 0, z = 0, in galactic cylindrical coordinates) [46], 
which features an important contribution at small wavelengths of starlight photons, because 
of the higher density of stars in the galactic bulge. For the '10x60' region we adopt the 
ISRF for R = 4 kpc, z = from [46], while for the '10—20' strips we use the ISRF computed 
for distances above the galactic plane, z = 5 kpc, in [47]. In the latter, the SL and IR 
intensity is significantly lower, as expected, but still sizable [47]. Table 1 summarizes the 
relevant parameters and fig. 1 illustrates the result. 

Finally, for the Dark Matter galactic profiles p(r) we consider three different mod- 
els determined by numerical simulations. Recent, state-of-the-art computations seem to 
converge towards the so called Einasto profile [48, 49] 



PEm[r) = p^exp 



a V V r.o 



a = 0.17. (27) 



The Navarro- Frenck- White profile [50] and the isothermal profile [51] 

PNFw(r) = Ps- (l + -] (28) PisoT(r) = , (29) 

represent instead previously standard choices. The values for the parameters and ps of 
the three models are given by 



DM halo model 


Ts in kpc 


ps in GeV/cm^ 


NFW [50] 


20 


0.26 


Einasto [48, 49] 


20 


0.06 


Isothermal [51] 


5 


1.16 



(note that all profiles are normalized at Pq at the location of the Earth). They are plotted 
in fig. 1. These profiles sensibly differ at the Galactic Center (from the cored isothermal 
profile to the more cuspy NFW). They are however also appreciably different a few kpc 
away from it. Since we consider integrated signals from regions of observation that extend 
along several kpc from the GC, a dependence on the choice of profile is therefore to be 
expected. 



4 Results 

In figure 2 we show some examples of the ICS signals (e^ ■ d^^ / deiAVL) for selected DM 
models, one for each of the observational regions that we study. 

In fig. 2a we plot the ICS gamma ray flux in the '5x30' region, from a DM candidate of 
^DM = 3 TeV annihilating into r+r" with (dann'^) = 1.7-10"^^ cm^/sec. We have assumed 
an Einasto DM profile. We plot the individual contributions of ICS signals onto the CMB 
(green), IR (red) and SL (blue) photon fields, and the resulting total flux (black thicker). 
As we see, the total flux skims through the EGRET data points at the largest energies. Our 
result compares with the result in ref. [24], obtained with GALPROP, modulo a slightly 
different choice for the Einasto profile and possibly the precise definition of the observation 
region. We find however a good agreement. 

In fig. 2b we plot the ICS gamma ray flux in the '10x60' region, from a DM candidate of 



10 



5x30 region 



10x60 region 



10-20 region 




Figure 2: Some examples of ICS signals from selected DM models, in the three different 
regions of the sky that we consider, superimposed to the relevant datasets. Left; in the 
region '5x30', EGRET datapoints and the signal from a 3 TeV DM candidate annihilating 
into T^T^ with (Jg^^^v = 1.7 ■ 10~^^ cm^/sec, choosing the Einasto DM profile of eq. (27) 
(in analogy with [24], where however the Einasto profile differs slightly). Center; in the 
region '10x60', EGRET datapoints and a signal from a 1.5 TeV DM candidate annihilating 
into with Uann'^ = 5 • 10"^'^ cm^/sec, choosing the isothermal DM profile of eq. (29) 

(in analogy with [26], where however the isothermal profile differs slightly). Right; in the 
region '10—20' strips, EGRET datapoints and preliminary FERMI datapoints, together 
with a signal from a 10 TeV DM candidate annihilating into W~^W~ with aannV = 5 ■ 
10"^^ cm^/sec, choosing the NEW DM profile of eq. (28). 



^DM = 1-5 TeV annihilating into with (cann'?^) = 5-10~^^cm^/sec. We have assumed 

an isothermal DM profile. Here as well the flux lies just below the EGRET data points. 
Our result here compares with the result in ref. [26], obtained with GALPROP, modulo a 
slightly different choice for the isothermal proflle. We again flnd however a good agreement 
between our computation and the fully numerical one. 

Finally, in flg. 2c we plot the ICS gamma ray flux in the '10—20' strips, from a DM 
candidate of Mdm = 10 TeV annihilating into W^W~ with (aann'w) = 5 ■ 10~^^ cm^/sec,^ 
assuming an NEW DM proflle. The flux here lies well below the EGRET and FERMI 
datapoints. 

The comparisons with the full GALPROP results also allows us to validate the ap- 
proximations of our semi-analytical treatment. In particular, the diffusion processes that 
we have neglected can redistribute somewhat the population of that originate the ICS 
signal. In the case of the isothermal proflle the effect is expected to be minor, as the DM 
density is essentially uniform over a broad central region of the galaxy. In the case of peaked 
proflles (Einasto, NEW) the can be moved from the denser central regions towards the 
outside. This would slightly reduce the signal from observational regions centered around 
the GC (such as the '5 x 30' and '10 x 60' regions) and instead increase it for regions 
outside the GC (such as the '10—20' strip). Comparing with the fully numerical results 
presented in [24, 26], we estimate at about a factor of 2 the systematic error introduced by 
our approximations ^. 

^These parameters are characteristic of the Minimal Dark Matter model, see e.g. [52]. 
''Note added: Ref. [29] also in fact finds that the diffusion-less treatment produces fluxes that are good 
within a factor of two over the whole photon energy range in the '10-20' region, see their fig. 4. 
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Figure 3: The regions favored by PAMELA (green bands), containing in particular the 
areas favored by PAMELA+ATIC (yellow areas), are compared with the bounds from ICS 
secondary radiation. The first column of panels refers to DM annihilations into e~^e~ , the 
second into fj^~^fi~ and the third into t^t^ ; the three rows assume respectively an Einasto, 
an NEW and an isothermal profile. In each panel, the bounds from EGRET data in the 
'5x30' region are plotted with a short dashed red line, those from EGRET data in the 
'lOx 60' region with a solid red line and those from EGRET data in the '10—20' strips with 
a dashed red line. The preliminary FERMI bounds in the '10—20' strips are plotted with 
a dashed blue line. The ICS bounds are computed adopting the simplifying assumption of 
neglecting the diffusion of the source e^, see text for the full discussion. 
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Figure 4: ^45 in the previous fig. 3, hut for W^W~ , bb and ti annihilation channels. Since 
a DM particle fitting the PAMELA data has to be multi- Te V, the green hands start at larger 
masses. There is no possibility to fit the ATIC data in these channels. 
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We are only including the ICS contribution in these plots and in our work, and not 
the prompt 7-rays from DM annihilations. The latter ones are model dependent and are 
expected to be subdominant in these low average DM density regions, especially for large 
DM masses. 

Scanning systematically a large range of DM masses Mdm (100 GeV - 20 TeV) and 
considering all the annihilation channels of eq. (25), we derive then constraints on the an- 
nihilation cross section by the following conservative prescription. For each observational 
region with the corresponding data points, fixed an annihilation channel, a DM distribu- 
tion profile and a DM mass, we impose that the ICS signal must not exceed any of the 
experimental data points by more than 3a. This determines a maximum annihilation cross 
section. Notice that this is the most conservative possible approach, as we do not assume 
anything on the background. An alternative procedure would be to take the EGRET and 
FERMI signals to be of astrophysical origin and then impose the requirement that the 
signal from Dark Matter, summed to such astrophysical background, does not exceed the 
data points. This would make our final constraints much stronger. Since we do not want, 
however, to exclude more portions of the parameter space that it is justified by current 
data, we stick to the first prescription. 

Figure 3 and 4 shows our results. We superimpose the ICS constraints to the regions 
identified by the PAMELA data (green bands, the different lines represent different positron 
and electron propagation models). We also indicate the subregions that allow to fit the 
ATIC peak as well (in yellow), see [7, 54] for the whole analysis. 

Focussing first on the purely leptonic annihilation channels (fig. 3), one sees that, 
choosing the Einasto profile, the ICS bounds exclude a large portion of the PAMELA 
regions. In general, the most constraining bound is the one from the EGRET data points 
in the '10x60' region, both because they are characterized by smaller error bars and because 
their spectrum falls rapidly at large energies (see fig. 2b). 

It is worth pointing out that the exclusion lines of the '10 — 20' strips increase more 
steeply at large masses with respect to the 5 x 30 and 10 x 60 region ones. This is due to 
the fact that the data in such regions go up to 80 GeV (see figs. 2a and 2b), whereas in 
the '10 - 20' strips they reach only 6 GeV for EGRET and 9 GeV for FERMI (see fig. 2c) 
and thus the fluxes are less constrained. One expects that, when FERMI will provide data 
in all the regions of the sky and cover a bigger range of energy, it will have the power to 
probe the parameter space down to much smaller DM annihilation cross sections. If the 
EGRET 'GeV excess' will be ruled out in the other regions as well, this actually means 
that it can be expected that the exclusion lines from FERMI will become very stringent 
at large masses. 

One also sees that the DM DM —>■ e^e~ channel is generically more constrained than 
fi'^H~ or T~^T~ , something which has to do with the fact that the original e~^e~ population 
from such channel is concentrated at larger energies, and therefore produces higher energy 
ICS photons, where the datapoints are more constraining. 

Assuming different DM profiles relaxes slightly the bounds. For an isothermal profile, 
some portion around Mdm = 1 TeV are reopened. We notice that the Einasto profile leads 
to the most stringent bounds: this is just because such a profile predicts a larger DM 
density at radii around several kpc (see fig. la), which is where the regions that we are 
considering are located. The fact that the density prescribed by the NEW profile becomes 
larger than Einasto very close to the Galactic Center is not relevant for these broad field 
integrated signals. We also recall that the bounds relative to the isothermal profile are not 
expected to be affected by our no-diffusion approximation, so they can be considered as 
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the most robust. 

Fig. 4 shows that the same quahtative features are present in the 'hadronic' annihilation 
(including W^+VT") channels. Notice that the ICS bounds are only a factor of (9(few) higher 
than the corresponding purely leptonic ones, but the PAMELA regions are less constrained, 
because they are located at smaller cross section values (especially for large masses) for 
these channels.^ 

We specify that we are always assuming unit boost factor and no other particle physics 
enhancement (such as the Sommerfeld or resonance enhancement). If such factors are 
present, then it simply means that the quantity constrained on the vertical axis of our 
panels in fig. 3 and 4 is the product B x av or S x av (where B denotes the astrophysical 
boost factor and S an enhancement): the PAMELA regions and the ICS contours simply 
rescale together, as the effect responsible for the increase would be most probably present 
in both cases (the local galactic halo from which the PAMELA positrons come and the halo 
regions probed by EGRET and FERMI in gamma rays) in the same way. An exception to 
this is the possibility that a single local DM clump is responsible for the PAMELA signal; 
in that case the ICS signals from other regions of the halo would not contain the effect. 
This possibility looks however rather unlikely [53]. A detailed analysis of the distribution 
of substructures in the halo, and of the DM velocity dispersion, would be necessary to go 
beyond these qualitative statements (see e.g. [14]). 

5 Conclusions 

We have studied the compatibility of the interpretation in terms of DM annihilations 
of the excesses of electrons and positrons, recently claimed by CR experiments such as 
PAMELA and ATIC, with the observation of diffuse gamma ray fluxes, as they would be 
produced by Inverse Compton scatterings on interstellar photons of such large populations 
of energetic in the galactic halo. We first presented a review of the basic formalism 
needed to compute the ICS signal in the context of DM annihilation, and then discussed 
its concrete application. We have performed the computation of the expected ICS fluxes 
semi-analytically (as opposed to running numerical galactic codes). We have worked under 
the simplifying assumption of neglecting the diffusion of the electrons away from their 
production spot in the galactic halo. Moreover, we have adopted simplifying modelizations 
of the Interstellar Radiation Field, based however on the precise computations in the 
literature. We find good agreement between the results of our procedure and the fully 
numerical computations. We can quantify at about a factor of 2 the systematic error 
introduced by our approximations. 

We have considered three main regions of observation and compared the expected fluxes 
with the existing EGRET data and with the preliminary FERMI points, deriving con- 
straints on the plane MoM/cann^^- The figures 3 and 4 summarize our results. They show 
that, for purely leptonic annihilation channels, large portions of the PAMELA regions at 
Mdm bigger than a few hundred GeV are excluded by a factor ranging from 'a few' up to 
about one order of magnitude, if an Einasto DM profile is assumed (for an NEW profile a 
very similar situation holds). This includes in particular the Mdm ~ 1 TeV region. If an 

^At small masses, the exclusion line from the preliminary FERMI datapoints in the '10—20' strip 
becomes less stringent than the corresponding line of EGRET. This is due to the fact that, for the 'hadronic' 
annihilation channels, the primary electron and positron spectra, and therefore the ICS 7-ray spectra, are 
broader, and in particular are sizable at low energy, where EGRET has currently more points than those 
available in the preliminary FERMI release that we are using. 
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isothermal profile (disfavored however by numerical simulations) is chosen, the constraints 
are somewhat relaxed, but still most of the Mdm > few TeV area is excluded. For the 
hadronic and W^W~ annihilation channels the qualitative results are similar, but larger 
portions of the PAMELA regions remain allowed. 

Complementary constraints on DM annihilation are imposed by high energy gamma 
rays (generated directly from the DM annihilation process) from the galactic center region 
and from satellite galaxies and by synchrotron radiation (generated by in the galactic 
center's magnetic field). A model-independent analysis, analogous to the one that we 
performed here, has been carried out in [54]. Its results show that the regions of the 
parameter space that allow to fit the PAMELA (and ATIC) data are disfavored by about 
one order of magnitude if a benchmark Einasto or NEW profile is assumed. But choosing a 
smoother profile and / or assuming that a part of the cross section is due to an astrophysical 
boost factor that would not be present in dwarf galaxies and the Galactic Center due to 
tidal disruption re-allows part of the space. 

The synchrotron and gamma ray constraints are therefore qualitatively quite similar and 
quantitatively somewhat stronger than those imposed by ICS, studied in this paper. They 
are however apparently quite dependent on the choices of the DM distribution, as expected 
from the fact that they come mainly from signals originating at the Galactic Center. On 
the contrary, the ICS signals (and therefore the constraints) are more robust with respect 
to these DM-distribution related issues. 

To conclude, our results show that the ICS gamma ray fluxes have the power of ex- 
cluding large portions of the parameter space (in DM mass and annihilation cross section) 
individuated by the PAMELA data, especially for masses larger than about 1 TeV, for 
purely leptonic annihilation channels and for benchmark Einasto or NEW profiles. This 
adds tension to the DM interpretation of the PAMELA data, on top of the already exist- 
ing constraints from prompt high energy gamma rays and synchrotron radiation from the 
same DM annihilations [54]: gamma rays seem to disfavor the DM interpretation of the 
PAMELA data in a multiwavelength range. Euture results from the EERMI satellite will 
have the power to probe larger portions of the parameter space. 
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